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Spatially Extended [P II] 1.188 fiui and [Fe II] 1.257 fiui Emission Lines 
in a Nearby Seyfert Galaxy NGC 1068 Observed with OAO/ISLE 
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Abstract 

We present J-band long-slit spectroscopic observation of NGC 1068 classified as a Seyfert 2 galaxy. 
J-band observations with OAO/ISLE provide clear detection of spatially extended [Fe II]1.257/im and 
[P II]1.188/im lines. We found that [Fe II]1.257/im/[P II]1.188/^m increases with distance from a central 
continuum peak. Observed line ratios around the nucleus (continuum peak) are consistent with a typical 
value expected from photoionization models, while the ratios at 3" — 4" (210 — 280 pc) east and west 
of the nucleus are slightly higher than this. In the off nucleus region of NGC 1068 we also found a 
possible association between [Fe II]1.257/zm/[P II]1.188/zm and the radio continuum. This suggests a mild 
contribution of shock ionization induced by a radio jet outside nucleus while photoionization by the central 
energy source is dominant near the nucleus. 
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1. Introduction 

The narrow line regions (NLRs), which extend to sev- 
eral hundred or kilo parsec scale around the galaxy cen- 
ter, are the exclusive structure of active galactic nuclei 
(AGN) where the spatially resolved observations are pos- 
sible. Therefore NLRs are often investigated as an impor- 
tant tool to study the ionization state of the interstellar 
medium (ISM) and/or chemical evolution in galactic scale 
(e.g., Nagao et al. 2006). Although it is widely accepted 
that the NLR is photoionized by ionizing photons radi- 
ated from a central engine, the possibility of shock ion- 
ization induced by a jet in off nucleus regions cannot be 
excluded since ionization photons decrease with distance 
from a nucleus (e.g., Fu & Stockton 2007). Thus, how 
much the shock contributes to the ionization of NLRs is 
very important to our understanding of AGN structure 
and in examining the utility of NLRs as a tool to investi- 
gate galactic-scale phenomena. 

Furthermore, the shock ionization of NLRs is getting 
a lot more attention. Recent dramatic progress of the- 
oretical simulations and observational studies of galaxy 
formation and evolution allow a quantitative comparison 
between both sides. In this context, a serious problem 
has arisen, i.e., theoretical simulations predict too many 
massive galaxies due to long-duration star formations in 
contrast to early-time quenching of star formation in ob- 
served massive galaxies (e.g., Croton et al. 2006; Bower et 
al. 2006). This problem cannot be solved even if a nega- 
tive feedback effect on star formation activity caused by 



supernovae is involved and so AGN feedback effect is con- 
sidered as a potential candidate of a solution: a massive 
galaxy likely has a supermassive black hole at its galaxy 
center, and inflow of ISM to a supermassive black hole 
invokes its AGN activity which releases vast gravitational 
potential energy to ISM resulting in suppression of star 
formation activity (e.g., Scannapieco et al. 2005; Sijacki 
et al. 2007). However, how the AGN activity transmits 
its energy to ISM remains a mystery. One possible phys- 
ical mechanism of the AGN feedback is shock ionization 
of ISM, i.e., the AGN activity inputs its energy to ISM 
through a shock heating induced by a jet. 

In previous studies of NLRs, line-ratio diagnostics to 
distinguish between shock ionization and photoionization 
have been examined. Optical diagnostics, however, can 
hardly discriminate between the two mechanisms, be- 
cause optical NLR spectra predicted by photoionization 
and shock ionization models are very similar to each 
other (Dopita & Sutherland 1995; Dopita & Sutherland 
1996; Allen et al. 2008). The near infrared line ratio of 
[Fe II]1.257/xm/[P II]1.188/im is one of the most power- 
ful indicators to discriminate photoionization and shock 
ionization. Both lines have similar critical densities and 
excitation temperatures, i.e., this line ratio is roughly pro- 
portional to the ratio of gas-phase abundance of iron and 
phosphorous. In contrast, iron is a well known refractory 
species and is strongly depleted in dust grains, whereas 
phosphorous is a non-refractory species. Photoionization 
alone (including H II regions and NLRs excited by ionizing 
photons from young stars and AGN central sources, re- 
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spectively) is relatively incapable of destroying the tough 
iron based grains, while these are easily sputtered by 
shocks. The [Fe II]1.257^m/[P II]1.188Aim ratio, there- 
fore, is high 20) in fast shock-excited regions and 
low 2) in normal photoionized regions (Oliva et al. 
2001). The actual ionization state of NLRs would be de- 
termined by the combination of photoionization and shock 
ionization, and so the observed line ratios are expected 
to change with the locations in NLRs ranging from [Fe 
II]1.257/Ltm/[P II]1.188/^m - 2 to - 20. 

NGC 1068 is one of the nearest AGNs, which has a com- 
pact radio jet around the nucleus and spatially extended 
radio lobe (Wilson & Ulvcstad 1983). This radio structure 
well coincides with a morphology of NLR (Capetti ct al. 
1997). Therefore, NGC 1068 is an ideal object with which 
to research the spatial distribution of [Fe II]1.257/im/[P 
II]1.188^m. In this paper, we adopt 14.4 Mpc as a dis- 
tance to NGC 1068 and 1" corresponds to 70 pc. 



ionized by AGN radiation because even ionizing photons 
from AGN central source can hardly destroy the tough 
iron based grains in NLR. 

However if shocks exist the grains arc easily destroyed 
and gas-phase iron increases. As a result shock ionized gas 
represents high [Fe II]1.257//m/[P II]1.188/im ratio. If we 
assume that iron based grains are completely sputtered by 
shocks, the ratio becomes ~ 50 for solar abundance. This 
is similar to that measured in supernova remnants (e.g., 
> 20 for LMC-N63A and LMC-N49 reported by Oliva et 
al. 2001). 

The actual observed line ratio in NLR of AGN is ex- 
pected to be between ^ 2 and ^ 20 since the ionization 
state would be determined by the combination of pho- 
toionization and shock ionization if these exist as men- 
tioned in section 1. 

3. Observations and data reduction 



2. Theoretical [Fe II] 1.257/xm/[P II]1.188/xm ra- 
tio 

Before describing the details of observation and data 
reduction, we summarize the theoretical background of 
[Fe II]1.257/im and [P II]1.188/im emission lines. 

The emission line intensity is proportional to the prod- 
uct of the density of the ion responsible for the emission 
line process (n^) and the electron density (ng), multiplied 
by a function / giving the rate of the process. Thus the 
intensity ratio of two emission lines radiated from ion 1 
and 2 is written as 



/(Al) _ 71i,l./l 



/(A2) ni,2/2 



(1) 



assuming the same spatial distribution of both ions 
(Osterbrock 1989). The function / involves the rate of 
emission line photons in the radiative transition from the 
excited level to the ground level written as 
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Here Aiq is the radiative transition probability, and qoi 
and qio are coUisional excitation and dcexcitation rate, 
which include collisional strength (f2). Thus, / can be 
calculated if A and ft are known. Oliva et al. 2001 derived 



n(Fc) < n{Fc+) 
n{P) n(P+) 



/([Fe II]1.257Aim) 
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using the collision strengths and transition probabili- 
ties of [Fe II]1.257/im and [P II]1.188/im (Krueger & 
Czyzak 1970; Zhang & Pradhan 1995; Mcndoza & Zeippen 
1982; Nussbaumer & Storey 1988) . For solar abundance 
ratio of ri,(Fe)/ri,(P) ^ 100 and typical depletion factor of 
Fe (-^ 0.01) and P (^ 1.0 because of refractory species), 
equation (3) gives [Fc II]1.257/xm/[P II]1.188/im close to 
unity, although the depletion factor of iron differs from 
object to object. Actually, [Fe II]1.257^m/[P II]1.188//m 
is < 2 in normal photoionized region, e.g., ~ 2 in Orion 
Bar (Walmsley et al. 2000). This is also true for NLR 



Long-slit spectroscopy was carried out from November 8 
to 12 2009 with ISLE (Yanagisawa et al. 2006; Yanagisawa 
et al. 2008), which is a near- infrared imager and spectro- 
graph for the Cassegrain focus of the 1.88 m telescope at 
Okayama Astrophysical Observatory (OAO). The camera 
used for the spectroscopic observations has a projected 
scale of 0".25/pixcl. The spectrum of NGC 1068 was ob- 
tained with a slit of 2".0 (= 8 pixels) width and J-band 
grating which yields a 1.11 — 1.32//m spectrum with a dis- 
persion of 0.166 fim/pix. The spectral resolution is ~ 
1300 measured from an OH emission line at the central 
wavelength. The slit was oriented to E-W (i.e., position 
angle = 90°) and centered on the J-band continuum peak 
of NGC 1068 (Fig. 1). We note that the position angle is 
fixed to 90° in ISLE spectroscopic mode. Therefore, the 
slit was not placed along the major axis of NLR and the 
direction of the radio structure, at a position angle of ~ 30 
° (Das et al. 2006; Crenshaw et al. 2010). It lies outside of 
the nominal bicone of NLR and away from the axis of ra- 
dio emission. Unfortunately other slit positions north and 
south of the nucleus could not be completed due to bad 
weather conditions. The acquisition consisted of a scries 
of two 2-minute exposures with the object set at different 
positions along the slit followed by dome flats and calibra- 
tion lamp of Argon and Xenon. The seeing size was 1".0 - 
2".0. Since the weather conditions were not good through- 
out the observation, we excluded poor data from a total 
6-hour exposure on source, resulting in the effective expo- 
sure time of 4.4 hours on source. The standard data reduc- 
tion was performed for all selected spectra, i.e., dark frame 
subtraction, flat fielding, wavelength calibration, and sky 
subtraction using IRAF software. To correct the atmo- 
spheric spectral response and the instrumental efficiency, 
spectra of NGC 1068 were divided by spectra of some A- 
type rationing stars (HIP5310, HIP10795, HIP14077, and 
HIP22774) with the same airmass. In this reduction, the 
black body and Pa/3 absorption features of rationing stars 
were removed by spectral fitting with a black body func- 
tion and Voigt profile, respectively. The assumed effective 
temperatures of rationing stars are 8270, 7500, 8200, and 
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9230 K for HIP5310 (A3V), HIP10795 (A7V), HIP14077 
(A5V), and HIP22774 (AlV), respectively. 

4. Results and discussion 

The obtained 2-D spectra of NGC 1068 are displayed 
in Fig. 2. We detected [Fe II]1.257Aim and [P II]1.188Aini 
lines as well as Pa/3 and [S IX]1.252/ini. The spatial ex- 
tent of [Fe II]1.257/im and [P II]1.188Aim are ~ 14" and 
~ 7". These values are clearly larger than the typical see- 
ing size of 1".0 — 2".0. Thus, wc concluded that spatially 
extended [Fe II]1.257/^ni and [P II]1.188/im were success- 
fully detected. 

The line fluxes were measured from spectral fitting anal- 
ysis with IRAF spec fit task (Kriss 1994), assuming single 
gauss and underlying linear functions for each emission 
line. We summarized relative line fluxes normalized by [P 
II]1.188/^m in Table 1, extracted from central 2".0 region 
and cast and west neighbor regions. The detailed spatial 
distribution of [Fe II]1.257/im/[P II]1.188/^m line ratios is 
shown in Fig. 3 (a). 

Oliva et al. 2001 reported that the line ratio in the 
central ~ 2" region of NGC 1068 is about 1.5. Since 
this value corresponds with the photoionization scheme as 
mentioned above, they concluded that in the central re- 
gion most iron is locked into grains and shock excitation is 
not the primary origin of [Fe II] line emission. This expla- 
nation is relatively straightforward because there would 
be a large number of ionizing photons near the nucleus, 
that is enough to dominate the ionization of surrounding 
gases. Our measurement of [Fe II]1.257/im/[P II]1.188^m 
■~ 1.3 in the central 2" region is consistent with this value. 

However, this argument may not be valid at off nucleus 
regions. Wc found that [Fe II]1.257/im/[P II]1.188^m 
increases with distance from a central continuum peak. 
While observed line ratios around the nucleus are consis- 
tent with a prediction by photoionization models, the ra- 
tios at 3" — 4" east and west of the nucleus (~ 560 pc) are 
slightly higher than the typical value of [Fe II]1.257/xm/[P 
II]1.188/^m in the photoionized region. Although there 
are only two research efforts devoted to the spatial distri- 
bution of [Fe II]1.257/xm/[P II]1.188/^m ratios of NLRs 
(NGC 4151 by Storchi-Bergmann et al. 2009 and Mrk 
1066 by Riffel et al. 2010), similar results were reported 
in both cases. Storchi-Bergmann et al. 2009 found that 
[Fe II]1.257/im/[P II]1.188//m ratios are higher at ~ 130 
pc away from the nucleus of NGC 4151 6) than that in 
its nucleus 2). They also pointed out a possible spatial 
correlation between [Fe II]1.257/Ltm/[P II]1.188/im and ra- 
dio continuum structure. They suggested that shocks in- 
duced by a radio jet release the Fe locked in grains and 
produce an enhancement of the [Fe II] emission at off nu- 
cleus regions. Similarly Riffel et al. 2010 found that Mrk 
1066 presents [Fe II]1.257/im/[P II]1.188Aim - 3 at most 
locations within ^ 470 pc from the nucleus, but in some 
regions close to the borders of the radio continuum struc- 
ture this ratio reaches values up to 9.5. They concluded 
that shocks seem to play a more important role in these 
regions. 
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Fig. 3 (b) shows VLA 4.86 GHz flux density as a func- 
tion of distance from the nucleus of NGC 1068 extracted 
from a same slit aperture as our near-infrared observa- 
tion with OAO/ISLE. The radio data was obtained from 
NRAO Science Data Archive^. At off nucleus region 
of NGC 1068 we flnd a possible association between [Fe 
II]1.257/zm/[P II]1.188/xm and the radio continuum like 
NGC 4151 and Mrk 1066. The higher ratios at off nucleus 
of NGC 1068 is likely attributed to a mild contribution 
of shock ionization to ionized gases. This may indicate 
that the interaction between the jet and ISM forms an 
expanding cocoon which induces the shock waves propa- 
gating perpendicularly in the direction of the jet axis (e.g., 
Scheuer 1974), while photoionization by central engine is 
dominant near the nucleus. 

5. Conclusion 

The line ratio [Fe II]1.257/^m/[P II]1.188Aim in the 
near-infrared wavelength range is a useful tool with 
which to examine the dust destruction by shocks. We 
investigated spatial distribution of this ratio in NLR of 
nearby Seyfert galaxy NGC 1068 with OAO/ISLE. [Fe 
II]1.257/xm/[P II]1.188/xm near the nucleus is close to 
unity consistent with a previous observation and with 
a ratio in a normal photoionized region. This indicates 
that photoionization by ionizing photons radiating from 
a central engine is dominant near the nucleus. We 
found that the ratio increases with the distance from the 
nucleus, and is slightly higher at 3" — 4" east and west of 
the nucleus than ratios typical of a photoionized region. 
We also found a possible spatial association between [Fe 
II]1.257/xm/[P II]1.188/^m and radio continuum around 
^ 560 pc from the nucleus. These flndings suggest a 
higher contribution of shock ionization induced by a radio 
jet at off nucleus. Except for NGC 1068, recently the 
spatial correlation between [Fe II]1.257/im/[P II]1.188/im 
and radio continuum over the several hundred parsec 
scale has been reported for NGC 4151 and Mrk 1066. 
Applying this kind of research to a number of other AGNs 
is the clue to revealing ongoing AGN feedback phenomena. 
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Fig. 1. J-band image of NGC 1068 obtained with OAO/ISLE in our observation. The long-sht position (P.A. = 90°) is shown by 
two solid lines. 
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Fig. 2. 2-D spectra in J band extracted from central ±15" region (a) and continuum-subtracted spectrum (b). 
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Fig. 3. [Fc II]1.257/im/[P II]1.188/im line ratio (top) and VLA 4.86 GHz flux density (bottom) as a function of distance from a 
continuum peak. Arrows in the top figure are lower limits calculated from 3 cr noise level around undetected [P II]f.l88/im. 



Table 1. Relative line fluxes normalized by [P fl]1.188/jm 



Line ID 


East 3".0 


Central 2".0 


West 3".0 


[P II]1.188Aini 


1.0±0.29 


1.0±0.02 


1.0±0.08 


[S IX]1.252/im 


1.05±0.02 


1.01±0.02 


0.73±0.07 


[Fe II]1.257/im 


1.79±0.02 


1.33±0.05 


1.63±0.07 


Pa^ 


2.35±0.04 


2.89±0.05 


2.55±0.08 



Spectra were extracted from central 2".0 region and east and west neighbor 3".0 regions. 
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